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Abstract 

In  this  article  it  is  demonstrated  that  the  compressive  response  of  uni- 
directionally  reinforced  composites  may  initiate  in  a  micro-buckling  mode 
and  subsequently  switch  to  a  micro-kinked  configuration.  The  foregoing 
possibility  derives  from  a  mechanics  model  which  considers  initial  hbers 
misalignments,  non-linear  shear  response  of  the  matrix,  shear  deformable 
fibers  and  stcKhastic  fiber  spacings.  The  latter  non-uniformity  in  fiber 
spacings  plays  a  majOT  role  in  the  generation  of  fiber  kiiUcs. 

Various  features  of  the  compressive  deformation  and  failure  process  are 
exhibited  by  computational  examples. 
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Introduction 


The  compressive  response  of  composites  has  been  studied  extensively 
during  the  past  three  decades  by  many  researchers.  Many  of  the  above  studies 
are  noted  in  comprehensive  listings  and  review  articles  by  Shuart  (1985), 
Camponeschi  (1991),  Guynn  et  al.  (1992)  and  Piggott  (1993),  and  will  not  be 
detailed  here.  Sutfioes  to  say  that  experimental  data  on  compressive  response 
exhibit  substantial  scatter,  which  is  partly  attributable  to  variability  of  the 
loading  mechanisms  and  partly  due  to  its  sensitivity  to  random  material 
flaws  and  misalignments.  On  the  other  hand,  many  analytical  and 
computational  models  predicted  compressive  strengths  in  excess  of 
experimental  values.  This  shortcoming  led  to  an  ongoing  effort  to  construct 
more  complicated  models,  which  incorporate  additional  material  and 
structural  parameters,  to  achieve  closer  correlation  between  predictions  and 
data. 


An  intriguing  aspect  of  the  compressive  response  of  fibrous  composites  is 
their  failure  by  kinking.  In  most  circumstances  p>ost-failure  inspections 
reveal  the  presence  of  titled  bands  of  broken  fibers,  separated  along  oblique 
straight  lines  from  the  remainder  of  the  test  sample.  Since  the  formation  of 
these  so-called  kink  bands  does  not  appear  accord  with  a  buckling  mode  of 
failure,  most  existing  models  address  the  compressive  failure  of  composites 
either  as  a  buckling  or  as  a  kinking  phenomenon,  to  the  exclusion  of  the 
other.  These  models  can  be  grouped  as  follows: 


(1)  Models  which  consider  buckling. 

These  include  the  work  of  Rosen  (1965),  which  seems  to  be  the  first  article 
on  compressive  failure  of  composites.  Considering  "shear-mode  buckling" 

that  model  predicted  a  failure  stress  Ocr  =  Gm  /(l  -  ,  where  Gm  is  the  shear 

modulus  of  the  matrix  and  ^fthe  fiber  volume  fraction.  That  prediction  is 
inadequate  for  two  reasons:  (a)  it  gives  Gcr  which  is  several  times  higher 

than  experimental  values,  (b)  the  relation  OcR~l/(l-%)  contradicts 
experimental  observations  which  show  that  Gcr  grows  linearly  with  (at 
least  up  to  -  0.55)  (e.g.  Piggott  and  Harris  (1980),  Morley  (1987)). 


Several  modifications  to  Rosen's  model  were  introduced  subsequently. - 

Primarily,  these  modifications  considered  non-linear  shear  response  of  the  - 

matrix  and  initial  fiber  waviness  (e.g.  Wang  (1978),  Lin  and  ^ang  (1992),  ^ 

Guynn  et  al.  (1992),  Highsmith  et  al.  (1992)  and  others  listed  in  the  ^ 

aforementioned  review  articles).  Additional  modifications  included  the  _ _ 

incorporation  of  fibers'  shear-deformation,  such  as  by  Davis  (1975),  or  the _ 

accoimting  for  large  deformations  of  the  fibers  by  Yin  (1992).  Though  the _ 

latter  model  stems  from  a  buckling  formulation,  it  is  worth  noting  that  it  •/ 
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proposes  a  criterion  for  kink  formation,  which  occurs  when  fibers'  curvature 
attains  a  critical  value. 

(2)  Models  whidi  consider  the  a-priori  existence  ctf  kinks. 

These  include  works  by  Evans  and  Adler  (1978),  Hahn  and  Williams 
(1986),  and  Budiansky  and  Fleck  (1992). 

In  all  the  above  works  the  tiber  reinforced  composites  were  viewed  as 
lamellar  regions  which  consist  of  Bber  and  matrix  layers  as  shown  in  Figure  1. 
It  should  be  noted  that  several  investigators  (Sadowsky  et  al.  (1967), 
Herrmann  et  al.  (1967),  Lanir  and  Fung  (1972)  and  Greszczuk  (1975)) 
considered  fibers  of  cylindrical  geometry.  All  the  latter  works  assumed  linear 
elastic  behavior  of  fibtf  and  matrix  materials. 

The  validity  of  representing  fiber-reinforced  materials  as  lamellar  regions 
was  questioned  recently  by  Weitsman  and  Chung  (1994),  where  severe 
disparities  were  shown  to  exist  between  strain  fields  in  lamellar  and  arrayed 
fibrous  geometries  subjected  to  buckling-mode  deformations.  This  issue  is 
the  subject  of  an  ongoing  investigation. 

The  main  purpose  of  the  present  article  is  to  present  a  model  for  the 
compressive  response  of  fiber  reinforced  composites  which  suggests  a 
transitional  mechanism  from  a  micro-buckling  form  of  deformation  to  a 
micro-kinking  mode  of  failure.  The  model  employs  the  lamellar  geometry  of 
Figiire  1  and  incorpoates  non-linear  shear  response  of  the  matrix  and  linearly 
elastic,  shear-deformable  fibers.  In  addition,  the  model  considers  two  kinds  of 
geometrical  imperfections,  initial  fiber  nusalignments  and  non-uniform  fiber 
spacing.  The  latter  consideration  is  the  novel  aspect  of  this  work.  In  closely 
spaced  fibrous  domains  the  narrow  matrix  regions  sustain  increased  levels  of 
shear  strains  which  approach  yield  or  failure  limits.  The  resulting  softening 
in  the  support  provided  by  the  matrix  to  the  fibers  is  shown  to  cause  localized 
fiber  failings,  accompained  by  overburdening  the  more  widely  spaced  fibers. 
Though  compressive  failure  still  occurs  by  buckling,  the  immediate  post- 
buckled  configuration  is  shown  to  consist  of  kinked  fibers.  The  sudden 
transition  from  buckling  to  kinking  is  due  to  the  abrupt  variation  in  the 
deformed  configuration  which  follows  the  buckling  instability.  This 
comprehensive  accounting  for  the  seemingly  disparate  phenomena  of 
budding  and  kinking  is  the  main  asset  of  the  present  model. 

Non-uniform  fiber  spadngs  were  considered  by  Chung  and  Weitsman  in  a 
previous  work  (1993).  However,  in  that  work  die  fibers  were  modelled  as 
Bemoulli-Euler  beams  and,  in  the  absotce  of  shear  deformations,  could  not 
develop  kinks.  It  is  interesting  to  note  that  under  the  constraints  inherent  in 
the  Bemoulli-Euler  beam  theory  the  absence  of  kinks  was  compensated  by  the 
presence  of  highly  concentrated  shear  forces  in  the  post-buckling  range  of 
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response.  An  abbreviated  version  of  the  current  article  with  computational 
results  based  upon  material  data  reported  by  Guynn  et  al.  (1992)  is  due  to 
a{^>ear  in  the  Proceedings  of  the  Twelfth  U.  S.  National  Congress  of  Applied 
Mechanics  (1.  Chung  and  Y.  J.  Weitsman,  1994). 


Basic  Equations 

Let  a  uni-directionally  reinforced  fibrous  composite  be  represented  by  a 
layered  medium  which  consists  of  hber  and  matrix  layers  of  thicknesses  2h 
and  2(c-h),  respectively  as  shown  in  Figure  1.  Let  ^  h/c  and  ^ni  ==  (c-h)/c 
denote  volume  fractions  where  here,  and  in  the  remainder  of  this  article, 
subscripts  f  and  m  refer  to  the  fiber  and  matrix,  respectively.  Let  x  and  y 
denote  coordinates  parallel  and  normal  to  the  layers,  with  corresponding 
displacements  u  and  v. 

Consider  an  initial  waviness  Vq  in  the  Hber  layer,  given  by 

VO  =  8oco6(icx/L)  (1) 


and  let  the  fibers  deflect  in  a  shear  mode  of  buckling  as  shown  in  Figure  2. 
(Rosen  (19d5),  Garg  et  al.  (1973)). 


The  length  L  prescribes  the  micro>buckling  length  associated  with 
compressive  loading  N  parallel  to  the  fibers.  We  assume  that  the  matrix  layer 
responds  in  shear  only,  sharing  a  common  displacement  v  with  the  fit^r, 
while  Ubi  varies  linearly  in  x  across  the  layer’s  thickness.  In  addition,  we 
suppose  that  the  deformation  within  the  fiber  layer  is  expressed  by 

Timoshenko's  shear-deformation  model,  where  ¥  denotes  the  independent 
rotation  of  the  planar  cross-section. 

Accounting  for  bending  efiects,  the  total  displacement  of  a  fiber  region  Uf 
is  given  by  (Washizu  (1975)) 


Uf 


Uf  +  ^^  [(v'  +  Vop-(vo'P]dx-yv 


(2) 


with  the  corresponding  strains 


(3a) 
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1<y«V'-W 


(3b) 


while 


•''’*2^)  [U.K2c-h).U,(x,hl) 


,  ^  dUf  ^ 

V  - ^  *v  +— y 

♦m  dy  ♦m 


In  equations  (2)  >  (4),  and  in  the  sequ^,  primes  denote  derivatives  with 
respect  to  X. 

Consider  linear  elastic  response  for  the  fibers 

=  Ef  Ex  /  ^  =  Gf  7^  (5) 

and  a  non>linear  shear  response  for  the  matrix,  scaled  by  the  initial  shear 
modulus  Gm 


^»Gs, 


The  total  relative  deflection  between  x  =  0  and  x  =  L/2  (say)  is  given  by 


4-u,(L/2)  +  l 


j  [('''  +  VoP-(''oP]' 


As  noted  in  the  Introduction,  random  fiber  spacing  accounts  for  an 
essentially  novel  feature  in  the  present  model.  Following  basic  statistical 
notions  concerning  the  subdivision  of  a  linear  segment  into  M  sub-segments 
of  random  lengths  (Dixon  and  Massey  (1969)),  we  assume  that  the  cell  sizes  2c 
are  distributed  according  to  a  Poisson's  point  process,  with  a  cumulative 
distribution  function 

P(C>c)®exp(-c/c)  (8) 
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In  equation  (8)  2c  is  the  nominal  average  of  the  cell  sizes.  Since  fiber 
regions  cannot  overlap,  namely  c  >  h  ("Gibbs  hard  core  process"),  equation 
(8)  is  modified  to  read 

P(C>c)-exp(-|^)  (9) 


The  cmresponding  probability  distribution  is 

P(c)..^«xp(-|^)  HO) 

In  the  sequel  we  shall  need  to  refer  to  the  probability  distribution  of  the 

fiber  volume  fractions  p  (^) ,  where  ^  ^  .  Employing  well  established  rules 
of  statistical  analysis  we  obtain 

pW-=^  -V  «pf-  M  <”) 

<hn  «  \  *f  I 


For  computational  purposes  the  distribution  function  p(^)  will  be 
discretized  and  expressed  by  a  finite  numbo*  of  Dirac  delta  functions. 


Field  Equationa  foi  Randomly  Spaced  Laycra  In  CompiBaaipn 

In  view  of  equation  (4),  and  subsequently  equation  (6),  it  follows  that  both 

iSy  and  are  functions  of  the  ranodm  variable  and  thereby  random  as 
well.  However,  to  ascertain  common  displacement  A  for  all  laya«  we  assume 
identical  displacements  u  and  v  for  all  layers  inspite  of  the  randomness  in 
their  widths  2c.  We  posit  that  the  effects  of  random  spacing  are  taken  up  by 

distinct  values  of  the  rotation  ¥,  which  vary  from  layer  to  layer. 
Consequently,  a  descretized  version  of  p(^)  which  corresponds  to  M 
values  of  (<n)  j  (i  «  1, . . .  M)  will  necessitate  the  consideration  of  M  values  of 
¥i  (i  - 1, . . .  M). 

The  above  assumption  involves  an  approximation  whose  validity  can  be 

estimated.  It  will  be  shown  below  that  a  solution  for  u,  v,  ¥i, . . .  ¥m  which  is 
based  upon  energy  considerations  for  the  entire  composite  fails  to  satisfy 
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equilibrium  conditions  for  any  of  the  M  individual  layers.  The  magnitudes 
of  the  imequilibrated  transverse  forces  (distributed  transverse  loads  qj(x),  i  »  1, 
.  .  .  M)  provide  the  required  estimate  of  error.  It  turns  out  that  the  above 
magnitudes  never  exce^  1%  of  all  calculated  stress  values,  attesting  to  the 
suitability  of  the  approximation. 


Application  of  the  principle  of  virtual  work  to  the  entire  composite 
(Washizu  (1975)),  upon  accounting  for  the  randomness  of  the  cell  sizes  c, 
yields. 


dc  +  N5A  *  0 
(c) 


(12) 


In  view  of  expressions  (3)  through  (6),  which  contain  the  random 

variables  and  and  since  by  hypothesis  ¥  depends  on  c,  it  follows  that  die 
expression  within  the  curly  brackets  on  the  left  side  of  equation  (12)  is  a 
function  of  the  random  variable  c. 


Consider  a  discrete  probability  distribution,  represented  by  Wf, 
commesurate  with  equations  (10)  and  (11),  which  corresponds  to  M  cells  of 
widths  2q  (i*  1, . . .  M).  In  this  case  the  integral  in  equation  (12)  is  replaced  by 
a  sum  and  we  obtain 


:  ?  Wi  f  j  {<^5e|[  +  tJy5)Jy)dy+  T ^ery^dy 

i-l  Jo 


dx  +  N5A  =  0 


(13) 


Note  that  and  contain  M  2  independent  variations, 

namely  5u,  5v '  and  5yi  (i  =  1, . . .  M). 

Substitution  of  equations  (2)  -  (7)  into  expression  (13),  integrating  by  parts 
and  collecting  terms  which  multiply  5u,  5v '  and  5vi  (i  =  1, . .  .  M)  we  obtain 
the  following  (M  +  1)  field  equations  for  v  and  ¥i  (i  =  1, . . .  M) 

-N(v’  +  V|;)  +  2h2  jGi(v'-Vi)  +(w^  GS,F[v'  +  (<)t/4jjvJ  j  =0  04) 
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- 2hGf(v'- vi)  +  2hGS,F[v  '  +  (<j>f /♦»), V J  =  0  (i  =  1, . . .  M) 


(15) 


with  the  boundary  conditions 

Vi(L/2)  =  0.  Vi(0)  =  0  (i=l....M).  aiidv'(0)  =  0.  (16) 

Note  diat  the  coupled  system  of  equations  (14)  and  (15)  is  non-linear  due  to 
the  presence  of  the  function  F,  whi^  is  non-linear  in  its  argument. 

In  equations  (15)  If «  h^/12.  Note  that  the  variation  5u  yields  the  trivial 
result  N  »  constant  »  Ngppijed/  require  no  further  consideration. 

Unlike  equations  (15),  which  apply  to  the  M  individual  cells,  equation 
(14)  expresses  the  lateral  equilibrium  of  the  entire  composite.  To  ascertain 
lateral  equilibrium  of  each  cell  it  is  necessary  to  impose  distributed  transverse 

loads  q|(x)  (isl, . . .  M)  and  consider  their  virtual  work  I  qjSvdx. 

Reemployment  of  the  principle  of  virtual  work  for  each  individual  cell, 

collecting  terms  that  multiply  6v  and  utilizing  the  solution  for  v  and  Vi 
which  corresponds  to  the  system  of  equations  (14)  -  (16),  we  get 

M 

qi  =n(v  *'  +  Vq')-  2h2  jCflv"'- vj)  +  (6ta/^GSiF'[v'  +(^/^Bi).V.]  j  (17) 
i»l 


The  relative  magnitudes  of  serve  as  measives  for  the  validity  of  the 
approximation  inherent  in  assuming  common  u  and  v  in  all  layers.  It  will  be 
shown  in  the  next  section  that  for  polymeric  composites  all  the  above  ratios 
appear  to  be  less  than 


A  non-dimensional  version  of  equations  (14)  and  (15)  is  obtained  upon 
introducing  X  =  x/L,  V  =  v/L  and  the  non-dimensional  factors  Of  =  , 

tiflf 


„2  _ 

E* 


EA  • 


Then  upon  denoting  Y  »  dV/dX,  equations  (14)  and  (15)  read 


8 


0  (18) 


<4  Vj  +  (♦m/^  oS,  F  [y  j » 


and 


>(i,'  +  o?(y-<|(J-aJ,F[Y  +  (*r/«4¥j-0  O-I,  •••«)  (19) 

with  the  boundary  conditions 

vj  (1/2)  =  0,  Vi  (0) «  0,  (i  =  1, . . .  M)  and  Y(0)  =  0.  (20) 


The  solution  to  equations  (18)  -  (20)  was  obtained  numerically  and  pertains 
to  graphite/PEEK  (APC-2)  composites.  The  interval  0  <  X  <  1/2  was 
subdivided  into  K  equal  sub-intervals  and  derivatives  were  expressed  by 
means  of  central  finite  differences.  The  discretized  system  of  equations  was 

solved  iteratively  to  obtain  values  of  V  and  ¥i  for  increasing  load  N  in  the 
pre-buckling  range.  In  accordance  with  observations  by  Johnson  et  al.  (1991) 

the  matrix  material  was  assumed  fail  in  those  regions  where  7^  ^  0.05.  The 
response  in  the  post-buckling  range  was  evaluated  by  prescribing  the  location 

X*  where  »  0.05  in  the  most  recently  failed  cells  and  evaluating  the 
corresponding  deflected  shape  V(X)  and  applied  load  N. 

Complete  details  are  given  in  the  Appendix. 

The  computations  employed  values  of  Ef  =  214  GPa  and  Gf  =  13.8  GPa, 
which  reflect  the  significant  transverse  isotropy  of  the  AS4  fibers  (Aboudi 
(1991)).  The  in-situ  shear  stress-strain  response  of  the  PEEK  resin  was  based 
upon  the  reduction  of  data  collected  by  Kyriakides  and  Liechti  (1993).  The  in- 
situ  data  are  shown  in  Figure  3,  where  they  are  contrasted  with  data  for 
unreinforced  PEEK. 


The  in-situ  shear  data  were  expressed  by  the  relation  ^  ^  where 

A  s  3096  MPa,  B  »  169.93  MPa,  q  =  0.23781  and  t  in  MPa.  As  noted  earlier, 
shear  failiure  was  assumed  to  occur  at  a  strain  of  Yu  -  0*05*  In  addition,  we  took 


9 


L  «  400d  and  Sq  «  4d,  as  suggested  by  Kyriakides  and  Liechti  (1993),  with  a  fiber 
diameter  d  »  7.6  ^m. 

The  computational  scheme  considered  four  random  cell  sizes,  namely 
four  discrete  random  values  of  which  accord  with  the  distribution 
function  of  equation  (11).  These  are  shown  in  Figure  4. 

Computaitonal  results  are  shown  in  Figures  5-12. 

Figure  5  exhibits  the  variation  of  the  applies  stress  o  vs.  the  non- 
dimensionalized  maximal  lateral  deflection  v(0)/h.  The  short  plateau-like 
region  at  1.5  <  v(0)/h  <  1.8  is  due  to  matrix  failure  which  initiates  and  spreads 
within  the  cell  with  highest  fiber  volume  fraction  =  0.975).  Failure  by 
buckling  occurs  at  o  >  1220  MPa  (point  A).  In  a  displacemnt  controlled 
experiment  a  further  increase  in  A  will  be  accompanied  by  a  drop  in  o  and  a 
corresponding  reduction  in  v(0)/h.  The  post  budding  response  is  represented 
by  the  segment  AB  in  Figure  5. 

Results  for  the  non-dimensionalized  displacement  V  and  slope  Y  are 
plotted  vtf  sus  the  non-dimensional  distance  X  in  Figures  6  and  7  for  various 
load  levels  in  both  pre  and  post  buckling  ranges.  As  may  be  expected,  these 
curves  are  approximately  proportional  to  cos  X  and  sin  X  in  the  pre-buckling 
range,  where  deformed  shapes  are  dominated  by  the  form  of  the  initi^ 
misalignment.  Note,  however,  the  substantial  deviations  in  the  deformed 
configurations  within  the  post-buckling  range,  where  they  no  longer 
resemble  the  cos  X  shape  of  the  initial  misalignment.  The  most  significant 
departiues  occur  near  X  =  0.5,  where  matrix  failures  initiate. 

The  distributions  of  versus  X  are  plotted  in  Figures  8(a)  through  8(d). 
These  distributions  are  shown  within  the  four  distinct  cells  and  for  various 
levels  of  applied  stress.  Figure  8(a)  corresponds  to  a  pre-buckled  stress  level 
o  >  654  MPa,  Figure  8(b)  corres{X>nds  to  the  buckling  stress  a  =  1220  MPa  (point 
A  in  Figure  5)  when  the  matrix  within  the  cell  with  =  0.975  failed  over  the 
range  0.16  <  X  <  0.5. 

Hgiures  8(c)  and  8(d)  correspond  to  post-buckling  loads  of  a  *  730  MPa  and 

o  s  722  MPa,  respectivdy.  The  latter  stress  corresponds  to  point  B  in  Figure  5, 
which  lies  directly  below  point  A  in  that  figiue.  Note  the  spreading  of  the 
failure  region  in  the  cell  with  =  0.8  between  the  post-buckling  stages 
depicted  in  Figures  8(c)  and  8(d).  Namely,  matrix  failure  which  was  confined 
to  0.49  <  X  <  0.5  at  o  =  730  MPa  expanded  over  0.45  <  X  <  0.5  as  the  p>ost 
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buckling  load  level  proceeded  to  drop  to  o  =  7.22  MPa.  There  were  also 
increases  in  within  the  unfailed  cells  with  »  0.45  and  »  0.175. 

'Hie  central  theme  of  the  present  article  is  demonstrated  in  Figures  9(a) 

through  9(d).  Those  figures  exhibit  the  fiber  shear  strains  versus  the  non- 
dimer^nal  distance  X  for  the  same  circumstances  as  in  Figures  8(a)  -  8(d). 
The  essential  feature  of  those  figures  develops  within  the  post-buckling 

range,  where  discontinuities  in  are  iu>ted  to  occur  simultaneously  within 
all  cells.  These  discontinuities  are  located  at  X  »  0.49  in  Figure  9(c)  and  at 
X  «  0.45  in  Figure  9(d).  No  such  discontinuities  exist  in  the  pre-buckling  stage 
or  at  buckling. 


We  suggest  that  the  foregoing  discontinuities  in  Hiiy  portend  the  onset  of 
fiber  kinking. 

Rgure  10(a)  exhibits  the  computed  values  of  the  lateral  loads  q|()0  versus 
X  within  all  four  cells  at  buckling^  namely  at  o  =  1220  MPa.  These  loads, 
whidi  were  evaluated  through  equation  (17),  are  much  smaller  than  Xu  >  70 
MPa,  ati^fing  to  the  validity  of  the  premises  of  the  present  model.  For 
purposes  of  comparison,  the  same  lateral  loads  are  plotted  versus  X  in  Figure 
10(b)  when  the  fibers  are  modelled  as  Bemoulli-Euler  beams,  employing  a 
previous  analysis  by  Chung  and  Weitsman  (1993). 

Additional  effects  of  random  fiber  spadngs  are  exhibited  in  Figures  11  and 
12.  Figures  11(a)  and  11(b)  concern  the  evolution  of  fiber  curvatures  with  load 
within  the  pre-buckling  range  and  demonstrate  the  contrast  between 
iiniformly  spaced  and  randomly  spaced  circumstances.  It  may  be  noted  that 
the  randomly  spaced  case  gives  rise  to  highly  concentrated  curvatures,  with 
magnitudes  which  exceed  by  fiuee  or  foiu*  folds  the  levels  which  correspond 
to  the  uniformly  spaced  case.  If  one  accepts  the  premise  of  Yin  (1992)  that 
kinks  occur  when  fiber's  curvatore  exceeds  a  certain  threshold  level,  then  the 
current  analysis  indicates  the  existence  of  an  enhanced  likelihood  of  kinking 
due  to  non-uniform  fiber  spacing.  Such  kinking  may  preceed  failure  by 
buckling. 

Finally,  the  efiect  of  random  fiber  spacing  on  compressive  strength  is 
shown  in  Hgure  12,  where  buckling  failure  loads  are  plotted  versus  the  fiber 

volume  fraction  The  results  for  the  random  case  were  computed  for 

average  values  ^  (with  the  distribution  shown  in  Figure  4)  that  are  identical 
with  the  constant  for  the  uniformly  spaced  case.  Note  that  random  fiber 
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spadngs  result  in  lower  strengths  with  ^  which  accords  with  experimental 
observations  by  Piggott  and  Harris  (1980). 


ConduaipM 

A  mechanics  model  was  presented  for  the  compressive  response  and 
failure  of  uni-directionally  reinforced  polymeric  composites  loaded  parallel  to 
the  fiber  direction.  The  analysis  account^  for  the  non-linear  shear  response 
of  the  resin,  including  its  ultimate  shear  strain,  and  for  shear  deformable 
fibers.  The  model  incorporated  two  kinds  of  geometric  imperfections, 
namely,  initial  fiber  waviness  and  random  fiber  spadngs.  Excepting  the 
works  by  Chung  and  Weitsman  (1993,  1994),  the  latter  k^d  of  imperfection 
has  not  been  considered  elsewhere. 

Computational  results  were  evaluated  within  both  pre  and  p)ost-buckling 
ranges  of  compressive  response.  The  paramount  result  of  this  work  is  an 
indication  that  a  band  of  discontinmty  in  the  fibers'  shear  strains  can  occur 
immediately  beyond  failure  by  buckling.  Such  a  discontinuity  may  portend 
the  formation  of  kink  bands  which  are  observed  in  many  failed  specimens. 

In  addition,  it  was  shown  that  an  accotmting  for  random  fiber  spadngs 
reduces  the  predicted  values  of  compressive  strength  and  results  in  a 
correlaiton  between  the  above  strength  and  fiber  volume  fractions  which 
concurs  with  experimentally  observed  trends. 
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Appendix 


dVi 

Upon  denoting  the  non-dimensionalized  governing  equations 


(18)  and  (19)  read 


-xV+Yo)  +  oJy  +  ^  w/-a?yi 

i-l  \ 


^  +  aS(Y-Yi)-oS,Fi  =  0 


+  a2,^iU0 


(A-1) 


dX 


(i  =  l,...M) 


The  first  equation  applies  to  the  entire  composite,  as  implied  by  the 
summation  over  all  the  Voronoi  cells,  while  the  second  and  third  apply  to 
the  M  individual  cells.  We  thus  have  2M-f  1  equations.  This  system  of 
nonlinear  coupled  ditferential  equations  can  be  converted  into  algebraic 
equations  by  means  of  finite  differences. 

Let  N  denote  the  number  of  nodes  along  the  buckling  span,  resulting  in 

N-1  intervals  of  length  AX  and  mid-nodes.  For  the  mid  node  between  n  and 
n-fl  th  nodes,  which  is  indicated  by  n-i-1/2  in  the  following,  the  finite 
difference  version  of  (A-1)  is 


-X^Yn+i/2+Yo^l/2)  +  a? Yn+l/2  +  2  W-a^^f^4,+l/2  +  oS,— F|^1/2)  =  0 

W  \  06  / 


AX 


+  aKYm^l/2  -  Vj^l/2)  -  CcS»F^n+i/2  =  0 


(A-2) 


AX 


=  ^1/2 


At  this  stage,  the  non-linearity  due  to  non-linear  matrix  response 
expressed  by  F(Yxy)  is  still  retained  wifiun  the  system.  This  non-linearity  can 
be  eliminated  by  considering  the  incremented  quantities  Y-t-5Y,  ^-*-54/  and 
y-i-5y  which  satisfy  equations  (A-2).  Linearization  is  then  achieved  upon 
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e)q>andizig  the  non-linear  temvs  in  Taylor  series,  and  retaining  the  first  order 
terms  oidy.  Then,  the  first  of  equations  (A-2)  results  in 


8Y11+1/2 


Ro+1/2  +  2  WiPtn+i/2&|^l/2 


M 


j  (A-3 


where 


M  /  Ami 

Rn+1/2  *  MYn+l/2+Yo,iM^l/2)  *  «? Yi»+i/2  +  X  H  “fVUM-1/2  ’  ^ - ^PIih^I/2 

W  \  Ab 


As  shown  in  equation  (A-3),  5Y,  which  is  common  to  all  cells,  can  be 
represented  explicitly  in  terms  of  ofiier  variables  Syi.  This  enables  the 
substitution  of  5Y  into  the  remaining  equations  in  (A-2).  The  substitution  of 
(A-3)  into  the  truncated  Taylor  expansions  of  the  last  two  of  equations  (A-2) 
yields  the  following: 


X  Yl^4H-l/2(5Vk,n+l'*^Vly»)"^^^^^(^MH-l+S\|rvn)“U^n+i/2  (A-4) 

AX  ici  ^ 


8^1+S^  .  Svu»+1-Svtn  „ 

- r - + - »Vi^l/2 

2  AX 


where 


S|^j/2  «  -  -  c^Yn+i/rVtn+1/2)  +  amFMH-i/2 

AX 
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Tjaii+1/2 


Wkl\iKl/2PuH^V2 

/M 

2  2  WiQi^l/2-X 


UgM.i/2  =  S|^i/2  - 


PUh-1/2R|h-1/2 

X  wjQu+i/2-X 
VM 


Vuh-1/2 


AX 


+  5111+1/2 


It  should  be  noted  that  equations  (A-4)  apply  to  each  individual  cell.  Thus, 

we  have  M  values  of  4i  Vi  node.  Collecting  54i  and  Syi's  and 

forming  the  vectors 


(Hl-(85i,8fe...85Mr 


|Sv)-{Svi,S<t^...SiKMr 


equation  (A-4)  can  be  expressed  in  matrix  notation  as 


(A-5) 


the  An,  Bn,  and  Cn  are  2M  by  2M  matrices,  and  Dn  is  a  2M  by  1  vector. 
Their  components  are  given  by 


B„(i,M+i)  =  - 


AX 


Bn(M+i,M-t-k)  =  Tk4^i/2 
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Bri(M+i,M+i)-Ty^l/2 

An(U)  =  -^ 


A4i,M+i)  = 


Qi4H-l/2 


AX 


A,^M+U)  =  - 


AX 


QU^1/2 


Ari(M+i,M+k)  =  Tit4^i/2 
An(M+i,M+i)  =  Ty^i/2  - 

c;<w)-^ 

C,(M+i,i)-  -L 

AX 

Eai)  =  Vu+l/2 
CaM+i)  =  Ui^l/2 


In  die  above  the  ranges  of  k  and  i  are  k^l, . . .  i-1,  i+l, . . .  M,  and  i=l, . . . 
M,  respectively. 

Denote  { {S^}  l^y}  T  (^l-  Then,  the  governing  equation  for  the  entire 
length  of  the  composite  is  obtained  by  collecting  equation  (A-5).  Therefore, 


[AJ  [Cj  0  • 

IMI 

(Dl) 

[Bd  [Ad  [Cd 
• 

i«j) 

1  • 

(P2) 

• 

• 

[Bn-J  [An-J  [Cn-z] 

1  {^-2} 

* 

(Dn-zI 

0  [Bn-i][An-i]. 

IDn-i) 

(A-6) 
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This  banded  algebraic  equation  can  be  readily  solved  by  means  of  the  LU 
decomposition  (Na  (1979)).  Since  the  solution  O  is  a  vector  of  increments,  it 
is  accumulated  at  each  ito'ative  step  to  compute  4  V-  At  the  end  of  each 
iteration,  the  norm  of  (b  is  calculated  to  check  for  convergence.  When  it 
becomes  sufficiently  small,  the  iteration  is  halted,  and  cumulative  values  of  % 

and  y  are  regarded  as  the  convergent  solution.  5Y  is  computed  by  using 
equation  (A-3)  and  accumulated  in  the  same  maxmer  to  obtain  Y.  Deflection, 
curvature,  stress  components  and  other  quantities  are  computed  by  means  of 
a  post  processing  subroutine. 

To  determine  the  buckling  load,  both  low  and  high  bounds  of  applied 
stress  and  ntimber  of  computation  into^als  are  inputed  in  the  program.  The 
highest  load  value  which  still  gives  a  convergent  solution  is  taken  to  be  the 
buckling  load. 

As  shown  in  Figure  5,  an  unloading  process  is  observed  in  the  post- 
buckling  regime.  In  contrast  with  the  presumed  linear  elastic  behavior  of 
fiber,  it  is  necessary  to  clearly  define  the  stiffness  of  the  matrix  during 
tmloading.  For  that  purpose,  isotropic  hardening  and  elastic  unloading  of 
matrix  were  assumed. 

The  post-buckling  analysis  employs  the  buckling  solution  as  a  starting 
stage,  which  establishes  all  the  initial  values  of  strains  and  displacements  that 
correspond  to  the  buckling  load.  Per*  The  post-buckling  computation 
postulates  a  location  X*  in  ffie  most  recently  failed  cell,  which  extends  the 
region  of  failed  matrix  in  that  cell  to  X*<X<l/2.  Consequently,  that  region, 
like  previously  failed  matrix  regions  imder  Per/  does  not  contribute  to  the 
shear  stiffness  of  the  composite.  Subsequently,  the  composite  is  subjected  to 
an  assiuned  load  level  P<Pcr  tbe  numeric^  scheme  is  employed  to  obtain 
values  of  all  corresponding  stresses  and  strains.  Note  that  wherever  the 
scheme  predicts  unloading,  it  is  programmed  to  follow  a  linear  elastic 

unloading  path  of  the  matrix.  The  value  of  y  at  the  foregoing  location  X*  is 
then  compared  with  Yu.  If  y(X*)  >  Yu  the  level  of  the  load  P  is  reduced 

iteratively  until  the  establishment  of  equality,  yiX*)  «  Yu  is  within  a  prescribed 
tolerance.  In  this  manner,  we  generate  the  load-deflection  curve  in  the  post 
buckling  regime. 
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Figure  TiUta 


1^.  1  The  fiber  reinforoed  composite  represented  by  a  layered  medium. 

Hg.  2  The  shear>mode  buckling  configuration. 

Fig.  3  Non-linear  shear  response  of  PEEK  at  21*0.  In-situ  response  is 
reduced  from  composite  AS4/PEEK  data  by  Kyriakides  and  Liechti 
(1993)  and  compared  with  the  behavior  of  neat  resin. 

Rg.  4  The  selected  discretization  of  the  random  distribution  expressed  in 
equation  (10)  into  four  distinct  values  of  with  average  value 

^  s  0.6. 

Fig.  5  Applied  stress  vs.  non-dimensional  maximal  lateral  deflection 
v(0)/h.  The  pre-buckling  stage  OA  is  followed  by  the  post-buckling 
stage  AB. 

Fig.  6  Plots  of  the  non-dimensional  deflection  V  vs.  the  non-dimensional 
distance  X  at  various  levels  of  applied  stress. 

Hg.  7  Plots  of  the  non-dimensional  slope  Y  vs.  the  non-dimensional 
distance  X  at  various  levels  of  applied  stress. 

Fig.  8  The  distribution  of  matrix  shear  stress  vs.  the  non-dimensional 
distance  X  within  the  four  distinct  cells  ( - =  0.175;  ■ 

♦f  =  0.45;  o  o  o  o  o  s  0.8;  ♦  ♦  ♦  ♦  ♦  =  0.975)  at  various  levels  of 

applied  stress.  8(a)  o  »  654  MPa  (pre-buckling),  8(b)  o  »  1220  MPa 

(buckling),  8(c)  o  »  730  MPa  (post-buckling),  and  8(d)  o  =  722  MPa 
(post-budding.  Point  B  in  Figure  5). 

Fig.  9  The  distribution  of  fiber  shear  strain  vs.  the  non-dimensional 
distance  X  in  the  foiu*  distinct  cells  (symbols  and  stress  levels  same 
as  in  Hgure  8). 

Rg.  10  Lateral  stresses  q{(x)  at  o  »  -  1220  MPa  vs.  the  non-dimensional 

distance  X  in  the  four  distinct  cells,  (a)  Shear  deformable  fibers, 
(b)  Fibers  deforming  as  Bemoulli-Euler  beams,  (o  -  0.25;  •  « 

0.45;  •••♦£  =  0.75;  - =  0.975). 
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Bg.  11  Non-dimensional  fiber  curvatures  (dY/dX)  vs.  non-dimensional 
distance  X  at  four  levels  of  applied  stress  in  the  pre-buckling  range 

( . 0-654  MPa;  - o  -  941  MPa;  •••••  o  -  1021  MPa; 

- o  -  1220  MPa  (buckling)),  (a)  Uniformly  spaced  fibers 

with  -  0.6.  (b)  Randomly  spaced  fibers  (acc<»rding  to  Figure  4)  with 

^-0.6. 


Fig.  12  Variation  of  compressive  strength  (-  buckling  load)  with  fiber 

volume  fraction  in  randomly  spaced  case).  Comparison 

between  uniformly  and  randomly  spaced  cases. 
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Fig.  1  The  fiber  reinforced  composite  represented  by  a  layered  medium. 


Fig.  2  The  shear-mode  buckling  configuration. 
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Fig.  3  Non-linear  shear  response  of  PEEK  at  21*C.  In-situ  response  is 
reduced  from  composite  AS4/PEEK  data  by  Kyriakides  and  Liechti  (1993)  and 
compared  with  the  behavior  of  neat  resin- 


0.1 75  0.45  0.8  0.975 


w, 

Fig.  4  The  selected  discretization  of  the  random  distribution  expressed  in 
equation  (10)  into  four  distinct  values  of  <|>f,  with  average  value  <|)f  =  0.6. 
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Fig.  5  Applied  stress  vs.  non-dimensional  maximal  lateral  deflection 
v(0)/h.  The  pre-buckling  stage  OA  is  followed  by  the  post-buckling  stage  AB. 
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Fig.  7  Plots  of  the  non-dimensional  slope  Y  vs.  the 
non-dimensional  distance  X  at  various  levels  of  applied  stress. 


Matrix  Shear  Stress  (MPa)  Matrix  Shear  Stress  (MPa) 


Fig.  8  The  distribution  of  matrix  shear  stress  vs.  the  non-dimensional 
distance  X  within  the  four  distinct  cells  ( - (|)f  =  0.175;  ■■■■■■■■  =  o, 

o  o  o  o  o  <|>f  =  0.8;  ♦  ♦  ♦  ♦  4  <|)£  =  0.975)  at  various  levels  of  applied  stress. 
8(a)  o  =  654  MPa  (pre-buckling),  8(b)  o  =  1220  MPa  (buckling), 

8(c)  o  =  730  MPa  (post-buckling),  and  8(d)  o  =  722  MPa 
(post-buckling.  Point  B  in  Figure  5). 
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Fig.  10  Lateral  stresses  qi(x)  at  o  =  =  1220  MPa  vs.  the  non-dimeiisional 

distance  X  in  the  four  distinct  cells,  (a)  Shear  deformable  fibers, 

(b)  Fibers  deforming  as  Bernoulli-Euler  beams. 

(o  <|»f  =  0.25;  •<!>;  =  0.45;  •••(!>£  =  0.75;  =  0.975). 
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Non-dimensional  fiber  curvatures  (dY/dX)  vs.  non-dimensional 
distance  X  at  four  levels  of  applied  stress  in  the  pre-buckling  range 

( .  0  =  654  MPa;  - o  =  941  MPa;  •••••  o  =  1021  MPa; 

- o  =  1220  MPa  (buckling)),  (a)  Uniformly  spaced  fibers  with  <l>f  =  0.6. 


(b)  Randomly  spaced  fibers  (according  to  Figure  4)  with  <>£  =  0.6. 


Fig.  12  Variation  of  compressive  strength  (=  buckling  load)  with  fiber 

volume  fraction  <|»/  ( =  <|>f  in  randomly  spaced  case).  Comparison  between 
uniformly  and  randomly  spaced  cases. 
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